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ABSTRACT 

The mass of super massive black holes at the center of galaxies is tightly correlated with the 
mass of the galaxy bulges which host them. This observed correlation implies a mechanism of 
joint growth, but the precise physical processes responsible are a matter of some debate. Here 
we report on the growth of black holes in 400 local galactic bulges which have experienced 
a strong burst of star formation in the past 600 Myr. The black holes in our sample have 
typical masses of 10 6,5 — 1O T,5 M and the active nuclei have bolometric luminosities of order 
10 42 — 10 44 erg/s. We combine stellar continuum indices with Ha luminosities to measure a 
decay timescale of ~ 3 00 Myr for the decline in star formation after a starburst. During the first 
600 Myr after a starburst, the black holes in our sample increase their mass by on-average 5% 
and the total mass of stars formed is about 10 3 times the total mass accreted onto the black 
hole. This ratio is similar to the ratio of stellar to black hole mass observed in present-day 
bulges. We find that the average rate of accretion of matter onto the black hole rises steeply 
roughly 250 Myr after the onset of the starburst. We show that our results are consistent with 
a simple model in which 0.5% of the mass lost by intermediate mass stars in the bulge is 
accreted by the black hole, but with a suppression in the efficiency of black hole growth at 
early times plausibly caused by supernova feedback, which is stronger at earlier times. We 
suggest this picture may be more generally applicable to black hole growth, and could help 
explain the strong correlation between bulge and black hole mass. 

Key words: galaxies: starburst, active, bulges 



1 INTRODUCTION 

The existence of a strong correlation between the masses of su- 
permassive black holes and the galactic bulges whe re they reside is 
strongly suggestive of a mechani s m of joint growt h ( Gebhardt et al.1 
l2000l : iFerrarese & Merrittll200Ql : iTremaine et alll2002h . However, 
the physical m echanism(s) responsible for this observed re lation 
remain unclear (IFerrarese & Fo rd 2005; Cattaneo et al.l l2009). The 
very different physical scales involved in building bulges through 
star formation and feeding a central black hole pose a great chal- 
lenge for the possible mechanisms of causal connection. Sev- 
eral theor etical models successfully reproduce the observed rela- 



tion (e.g. Isilk & Reesl[l9 98: Haeh nelt et aD 1 19981; [Hopkins et all 
2006: ICiotti& Ostriker 2007; Khalatyan et al. 200|), most invok 



ing some form of "feedback" by which the energy radiated during 
the phase of rapid accretion onto the black hole prevents further 
accretion and shuts down star formation. However, observationally 
there is little evidence for effective direct "mechanical" feedback 
caused exclusively by energy released during accretion onto the 
black hole. 

Clues to the nature of the relationship between the growth of 
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bulges and black holes have been provided over the last decade 
through the significant observational progress that has been made 
in understanding the demographics of the AGN population and 
its relationship to the properties of the host galaxies. In the lo- 
cal Uni verse, the Sloan Digital Sky Survey (SPSS) spectroscopic 
survey <York et al.l2000l : IStrauss et al.ll200llAbazaiian et al.l2009h 
has provided the statistics necessary to calculate volume averaged 
quantities and to reveal which populations of AGN and galaxies 
are the sites of t he ma jority of present day black hole growth. 
Kauff mann etaD <2003l) showed that significant growth of black 
holes is presently happening in galaxies with the masses and struc- 
tural properties of early-type galaxies, but with significant am ounts 
of recent or on-going star-formation. iHeckman et al.l J2004I) then 
found that the majority of low-redshift black hole growth occurs 
in high growth rate phases (L/LEdd greater than a few percent) 
which must be of short duration relative to the Hubble time. This 
is consistent with the fact that integrating the luminous output 
of quasars (high L/LEdd accretors) over all cosmic times gives 
a black h ole mass density sim ilar to that observed in the local 
Universe {Yu & Tre maine 2002). Comparing the integrated ongo- 
ing starformation to t he integrated ongoing black hole accretion, 
Heckman e7aD <2004l) found a value close to that expected from 
the Mbh — Mbuige relation, implying that bulge formation and BH 
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growth are still tightly coupled even today. However, rapid growth 
of the population of bulges and their black holes is occurring today 
only in systems of relatively low-mass (so-called cosmic downsiz- 
ing). 

One long-standing idea linking the evolution of bulges and 
black holes is that there is a strong evolutionary connection be- 
tween intense c entral burs t s of st ar formation and the growth of 
black holes. In IWild et al] |2007) we found that in fact the ma- 
jority of the most rapidly growing black holes at the present day 
exist in star-forming galaxies with bulges, but with relatively ordi- 
nary recent star formation histories and no indication of a recent 
strong burst of star formation. However, we also found that galax- 
ies which are undergoing or have recently undergone the strongest 
bursts of star formation in the sample showed an increased prob- 
ability of hosting a rapidly growing black hole. This link between 
strong starburst and black hole growth may have been more impor- 
tant at higher redshifts where gas was more abundant (potentially 
leading to larger and more frequent bursts of star formation), and 
where the majority of the black hole mass in the universe was ac- 
creted. 

Progress on understanding the link between star formation and 
black hole growth has also been provided by detailed observations 
of very nearby AGN which can probe spatial scales of 10s to 100s 
of parsecs (an order of magnitude smaller than the typical scale 
probed by the 3" SDSS fibre aperture). Evidence for young nuclear 
stella r populations is f ound in around half of local AGN studied 
(see lDavies et al.ll2007l for a recent summary of the literature). In 
the future such an approach may be extended to provide a system- 
atic investigation covering several orders of magnitude in bolomet- 
ric luminosity and black hole mass, and including starbursts with 
no current black hole growth. Of direct importance to this paper is 
the discovery of a possible delay between the peak of the starburst 
and AGN activity for a handful of nearby AGN ( Tadhunte r et al.1 
1 19961 : lEmonts et alJbOOd lGotoll2006l : bavies et alJl2007h . Such a 
delay may also be rela ted to the very low a ccretion rate observed 
in the Galactic center feaganoff et al .1120031) . despi te the presence 
of considerable wind ma t erial f rom young OB stars dPaumard et al.1 
2006). ISchawinski et al.l J2007h used the Lick indices and photom- 
etry of a morphologically selected sample of elliptical galaxies to 
reveal that those with emission lines dominated by star formation 
typically have more recent star formation than those with emission 
lines dominated by AGN. They interpreted this as a delay between 
starformation and AGN activity. However, they neglect the fact that 
more recent star formation will hide increasingly powerful AGN, 
which will significantly affect their result. 

In this paper we bridge the gap between the gross trends of 
stellar populations with AGN properties uncovered by SDSS, and 
detailed studies of the recent star formation in the nuclei of local 
AGN. We use the SDSS to select the galaxies with bulges which 
have recently undergone the strongest bursts of starformation in 
the local Universe. Our sample is complete for ages up to 600 Myr 
after the onset of the starburst, considerably longer than the typical 
timescale for black hole accretion events. From this sample, we 
construct a time-averaged view of black hole accretion subsequent 
to a star formation episode. 

In Section [2] we present the sample and the method used to 
measure the age of the starburst from the stellar continuum. In Sec- 
tion [3] we derive black hole accretion rates and instantaneous star 
formation rates from the emission line fluxes, and quantify the com- 
pleteness of our samples. The results are presented in Section|4]and 
discussed in Section \5\ Throughout the paper, we calculate lumi- 
nosities using Qm = 0.3, Qa = 0.7 and Ho = 70km/s/Mpc. 



log N 




4000A Break Strength 

Figure 1. The distribution of 4000 A break strength and B aimer absorption 
line strength for 70,000 bulge-galaxies (grayscale, in logarithmic number). 
The spectral indices are generated using a Principal Component Analysis 
of a model spectral library as described in Wild et al. (2007). The star for- 
mation histories of galaxies that lie in different regions of the diagram are 
described in Section|2] The starburst sample is overplotted as coloured dots, 
where colour indicates the age of the starburst. The dotted box indicates the 
region from which the starforming comparison sample is drawn. 

For the reader who is less interested in the technical aspects of this 
work we recommend reading the introductions to Sections 2 and 3 
(i.e. skipping the subsections), before moving onto the results and 
conclusions. 



2 A LOCAL STARBURST SAMPLE 

AGN (the observational manifestation of growing black holes) are 
found in galaxies with bulges, from early type spirals through to 
massive ellipticals (Ho et al. 1997; Dunlop et al. 2003; Kauffmann 
et al. 2003; but see Filippenko & Ho 2003 for an exception) . 
Therefore we select as our starting point a base sample of ^70,000 
galaxies from the Slo an Digital Sky Survey (SDSS) data release 7 
dAbazaiian et al .112009b . selected to have high stellar surface mass 
densities characteristic of bulges (Kauffmann et al. 2003). We re- 
strict the sample to galaxies with redshifts 0.01 < z < 0.07. At 
these redshifts the 3 arcsec aperture of the SDSS optical fibres sam- 
ples the inner <2kpc radius of the galaxies , so that ligh t from the 
stars in the bulge dominates the spectrum (Gadotti 2009). Broad- 
line (i.e. Type 1 or unobscured) AGN are not included in our sam- 
ple, as light from the central nucleus would contaminate our study 
of the stellar population in the host galaxy. Full details of the sam- 
ple selection criteria are given in the following subsection, in the 
remainder of this section we introduce the concept of our method- 
ology. 

In the integrated spectrum of a galaxy the different signatures 
of stars of different ages can be used to obtain information about 
a galaxy's recent star formation history. In Fig. Q] we begin by in- 
specting the distribution of two common stellar continuum indices 
for the 70,000 bulge-galaxies, the 40 00 A break s trength and the ex- 
cess B aimer absorption line strength dWild et al.l20070 . The major- 
ity of bulge-galaxies show no evidence of recent or current star for- 
mation, they form the "red- sequence" which lies on the right, with 
strong 4000A break strength indicating that the galaxies contain 
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predominantly old stars. There is a significant minority of bulge- 
galaxies that are forming stars, these have younger mean stellar 
ages and therefore weaker 4000A break strengths and form the 
"blue- sequence". A small number of bulge-galaxies are undergo- 
ing a "starburst" i.e. there has been a sharp increase in the galaxy's 
star formation rate over a timescale that is short (~ 10 8 years) 
in comparison to the age of the galaxy. They differ from blue- 
sequence galaxies only in the strength and/or short timescales of 
the star formation episode. These galaxies are identified by their 
unusually weak B aimer absorption lines, strong UV-blue continua, 
and weak 4000 A breaks i.e. spectra dominated by light from O/B 
stars. These objects lie in the lower left of Fig. Q] As the star- 
burst ages, the Balmer absorption lines increase in strength as 
the galaxy passes into the pos t- starburst phase foressler & Gunnl 
1 19831 : ICouch & Sharpies! 1 19871) . i.e. A/F star light dominates the 
integrated galaxy spectrum. These objects lie in the upper left 
of the figure. The shape of the left hand side of the distribution 
in Fig. Q] therefore describes the evolutionary track of a starburst 
galaxy, with time since the starburst increasing from bottom to top, 
and burst strength increasing from right to lef t dCouch & Sharpies! 
1987; Kauffmann et al. 2003; Wild et al. 2007). Thus the evolution 
of the strongest starbursts that occur in bulges in the local Universe 
defines the outer lefthand edge of the distribution. 

The unique advantage of our sample for timing the starburst- 
AGN connection is that the time since the onset of the starburst 
can be measured accurately using population synthesis models. For 
ordinary blue-sequence galaxies, the time since a recent, smaller 
change in star formation rate is essentially unconstrained using in- 
tegrated spectra alone, making a similar measurement of the rela- 
tive timing of star formation and black hole growth impossible with 
the data in hand. 

To sample the properties of AGN in different evolutionary 
stages of the starburst, we select a constant number of the strongest 
starburst galaxies per unit age. We cannot use our indices to mea- 
sure the age of starbursts to arbitrarily old ages, firstly because the 
light from weaker starbursts fades rapidly until they become in- 
distinguishable from ordinary starforming galaxies, and secondly 
because even the stronger starburst galaxies suffer from a "burst 
age-mass degeneracy" at late times. We therefore trade maximum 
burst age against sample statistics in order to select a complete sam- 
ple of 400 bulge-galaxies which have undergone a starburst in the 
last 0.6 Gyr. The starburst sample is shown as coloured dots in Fig. 

m 

In the following subsections we give further details of the sam- 
ple selection and the calculation of the age of the starburst. 

2.1 Details of the sample selection 

Our bulge-galaxy sample contains galaxies with stellar surface 
mass density /x* = 0.5M*/tt x r 2 > 3 x 10 8 M o /kpc 2 , where 
stellar mass (M*) is measured from the 5-band SDSS photometry 
(J. Brinchmann, http://www.mpa-garching.mpg.de/SDSS ) and ra- 
dius (r) is calculated from the z- band Petrosian half light radius. 
Although this selection does not necessarily isolate solely those 
galaxies with large bulge-to-disk mass ratios, this mass density cut 
is based on previous results that find that AGN in SDSS are pre- 
domina ntly located in host galaxies with high stellar surface mass 
density fcaufftnann et al.ll2003l : lHeckman et all2004h . 

We note that our stellar surface mass density cut excludes 
most, although not all, pseudo-bulges. Recently, Gadotti & Kauff- 
mann (2009) have shown that pseudo-bulges contain 4% of the 
black hole mass in the local Universe. Empirically, we find that 



SDSS galaxies with pseudo-bulges have a slightly lower limit on 
their stellar su rface mass density of /x* > 1 x lO 8 M /kpc 2 
taadottil2009h and therefore we have verified that our conclusions 
remain unchanged when we use this lower stellar surface mass den- 
sity limit. 

We impose a per-pixel signal-to-noise ratio (SNR) limit on the 
spectra of 8 in the g-band to ensure accurate measurement of the 
stellar continuum and emission lines. At these low redshifts and 
relatively high stellar masses of bulge-galaxies, only 4% of spectra 
have SNR below this cut and decreasing the limit from 8 to 6 does 
not alter our conclusions. 

We remove galaxies with stellar velocity dispersions 
<70km/s, below which the SDSS spectra have insufficient spec- 
tral resolution to measure accurate stellar velocity dispersions and 
therefore reliable black hole masses. This effectively places a lower 
limit on black hole mass of 1O 6 3 M (Tremaine et al. 2002), and 
we emphasise that our results directly pertain only to black holes 
with masses above this limit. As with any hard cut on sample prop- 
erties, this cut on stellar velocity dispersion may introduce a small 
bias into our results: new systems will enter the sample at later 
times as the black holes grow and cross the sample selection thresh- 
old. However, with typical black hole growth factors of 10% (see 
Section O the effect should be small. 

Finally, we remove a very small number (five) of candidate 
"dusty starbursts" which have very large Balmer decrements (>2.5 
times the intrinsic ratio, see below) and are identified as having 
Ha luminosities that are too high for their starburst age. Due to 
the obscuration of the youngest and hottest stars behind optically 
think dust clouds, these dusty but actively star forming galaxies 
may potentially have continuum spectral in dices which place them , 
incorrectly into the strong- starburst track (Poggianti&Wul bOOOL 
Wild et al. in prep). The inclusion or not of these 5 objects does not 
alter the conclusions of this paper. 

2.2 Measuring starburst age 

In order to boost the signal-to-noise ratio of traditional indicators of 
recent star formation history, in particular the measurement of the 
Balmer absorption lines, we use indices defined from a Principal 
Component Analysis (PCA) of the 4000A break region of the spec- 
tra. The full method is described in detail in Wild et al. (2007). To 
measure the age of the starbursts, we build a model library of 10 7 
star formation histories composed of an old bulge population with 
super posed exponentially decaying starbursts dBruzual & Chariot! 
I2003L Chariot & Bruzual, in prep). This star formation prescrip- 
tion, in addition to being intuitive for the sample at hand, provides 
a model library that fully covers the observed distribution of index 
values and naturally explains the observed shape of the distribution 
in Fig.Q] 

The bulge population is formed from a tophat starburst of 
0.5 Gyr width, evolved passively for 13 Gyr. The star formation his- 
tory of the model bulge population is unimportant, but the age is set 
to give a precise match to the position of the red- sequence in our 
stellar continuum indices. The metallicity of all models in the li- 
brary is fixed to 0.5 Z to agree with the position of the peak of the 
red-sequence for our galaxy sample. We note that this is slightly 
lower than generally measured for samples of SDSS galaxie s with 
masses > 10 10 M Q dPanter et al.ll2008l : iGallazzi et alJl2006k Each 
method uses different stellar continuum wavelength ranges and fea- 
tures, and thus the discrepancy may reflect a real difference in the 
metallicity of stars contributing most light to the bluer wavelengths, 
or a small inconsistency in the stellar population synthesis mod- 
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elling. Bursts of star formation with exponentially declining star 
formation rates (SFRoc exp^ _t//r ^) are superposed on this old stel- 
lar population, allowing burst mass fraction, burst age, decay rate 
and dust to vary with uniform priors. Burst age is allowed to range 
from to 1.5 Gyr, burst mass fraction from to 50% and dust (tv) 
from to 2. While the decline in star formation subsequent to a 
starburst is not constrained from the integrated spectrum of an indi- 
vidual galaxy, the entire distribution of continuum spectral indices 
observed in galaxies does provide some additional constraint. For 
a model in which star formation rate declines exponentially, r val- 
ues much greater than 0.3 Gyr are excluded as they wo uld not pro- 
duce the observed population of post- starburst galaxies ( Wild et al. 
2009). Very short values of r are equally excluded, because a burst 
that is strong enough to reproduce the starburst galaxies with short 
r, evolves along a track which lies to the left of the observed edge 
of the distribution in Fig.Q] where no galaxies are observed. Requir- 
ing consistency between the decay in Ha luminosities, which probe 
instantaneous star formation rates, and the age estimated from the 
continuum provides an additional constraint. Finally, we allow ex- 
ponentially declining starbursts with 200 < r/Myr < 350. We 
note that this starburst timescale is considerably longer than the "in- 
stantaneous" starbursts often assumed in the literature when mea- 
suring the ages of stellar populations in starburst and post- starburst 
galaxies. Our longer decay timescales result in slightly older de- 
rived ages, and may explain the population of galaxies with both 
strong B aimer absorption and emission lines (so-called e(a) galax- 
ies)Q. 

For each galaxy, the age of the starburst is estimated from 
the mode of the probability distribution function of model ages in 
the standard way. Starburst ages are well constrained, distributions 
are unimodal until beyond the maximum age presented in this pa- 
per, and ages are not sensitive to priors or dust reddening effects 
of the stellar continuum. Precise starburst mass fractions for indi- 
vidual objects are less well constrained than starburst ages. Rather 
than relying on estimated burst mass fractions to define our sam- 
ple of the strongest starburst galaxies, we select a constant number 
of galaxies with the lowest 4000A break strength at each starburst 
age. As described above, the total number of starbursts selected 
for the strong- starburst track sample trades completeness at late 
ages against sample statistics. At late ages, old/strong starbursts 
become degenerate with younger/weaker starbursts. Selecting 20 
starbursts per 30 Myr time bins gives a statistically adequate total 
sample of 400 galaxies complete to 600 Myr without risking the in- 
clusion of galaxies at later times which are in fact younger, weaker 
starbursts. 600 Myr is long compared to the expected timescale of 
AGN accretion events (10 7 to 10 8 years), and our results present 
"volume- averaged" quantities, averaging over a large enough sam- 
ple of galaxies to smooth out the stochastic nature of AGN activity. 



Figure 2. Trends of key galaxy properties with starburst age. In all pan- 
els, blue points indicate star-forming galaxies, orange indicate starforming- 
AGN composites, red indicate pure AGN and black open circles are unclas- 
sified due to having too weak lines. Typical errors are noted in the upper 
left of each panel where significant. From top to bottom as a function of 
time since the onset of the starburst: (i) total galaxy stellar masses, calcu- 
lated from the photometry and mass-to-light ratio taken from the likelihood 
distribution of the starburst models described in Section |Z2| (ii) black hole 
masses, measured from stellar velocity dispersion and assuming the M-cr 
relation ( Tremaine et al. 2002); (iii) stellar mass surface density, measured 
from the total stellar mass and the Petrosian half light radius in the z-band; 
(iv) concentration index, measured from the r-band Petrosian radii which 
enclose 50 and 90% of the light; (v) physical size, measured from the r- 
band Petrosian radius which encloses 90% of the light; (vi) redshift. 



2.3 Summary of sample properties 

To summarise our selection criteria: 

• Our base catalog is drawn from the SDSS DR7 spectroscopic 
catalog of spectroscopically identified galaxies with extinction cor- 
rected r-band Petrosian magnitude 14.5 < r < 17.7. 



More complicated star formation histories, such as those with an initial 
slow decline and fast later decline, could also produce the observed popula- 
tion of post-starburst galaxies. However, this would only result in an small 
overall stretching/contracting of the starburst age axis in places. 
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• We isolate galaxies with bulges based on their stellar surface 
mass densities ji* > 3 x lO 8 M /kpc 2 . 

• Our sample is limited to 0.01 < z < 0.07 to ensure the spec- 
trum observed through the 3" fibre aperture is dominated by the 
light from the bulge. 

• Due to the velocity resolution of SDSS spectra, we only accept 
galaxies with stellar velocity dispersion, measured within the fibre, 
of > 70km/s. 

• From this base catalog, we select the 400 bulges undergoing 
the strongest central starbursts, with starburst age <600Myr. This 
sample is complete in the sense that we select an equal number of 
starbursts per unit age. 

The starburst sample is shown as coloured dots in Fig. [TJ the 
smoothness of the line traced by its right-hand edge is impressive 
and supports our interpretation of the distribution of the spectral 
indices. The stellar masses, black hole masses, stellar mass den- 
sities, concentration indices, sizes and redshifts of galaxies in our 
starburst sample are shown in Fig. [2] as a function of starburst age. 
In Appendix [A] we present two panels of images of the youngest 
and oldest starbursts in our sample. Note in particular the size of 
the SDSS fibre through which the spectrum is taken indicated in 
the top left of the Fig.lATI 

The total stellar masses of the galaxies presented in Fig. [2 
range typically between 5 x 10 9 and 5 x 1O 1O M . These masses 
have been estimated from the probability distribution function of 
mass-to-light ratios output from the models described above, and 
are consistent with those determined independently from SDSS fi- 
bre magnitudes. Based on these model fits, our sample of starburst 
galaxies have converted a mass of gas into stars that is typically 
around 10-15 percent of the total stellar mass of stars in the fiber 
aperture. Because starburst mass fractions derived from stellar con- 
tinua are slightly prior dependent, we also estimate burst masses 
from the Ha luminosities of the galaxies (presented in the follow- 
ing section). These suggest a slightly higher typical burst mass frac- 
tion of 25 percent. 



3 BLACK HOLE ACCRETION RATE 

We are interested in two quantities that describe the growth of 
the black hole. Firstly, the black hole accretion rate (BHAR, 
dM a cc/dt) that we estimate directly from the extinction cor- 
rected [O III] line luminosity using a bolometric correction of 
600, which corresponds to the conversion BHAR/(M0/yr) = 

4 x 10 _10 Loin/L©. Secondly, the black hole growth rate 
(dMacc/dt/MsH) which characterises the timescale for the black 
hole to double its mass, and which we quote in units of the Ed- 
dington Luminosit)Q (LEdd), assuming a radiative efficiency of ac- 
cretion of 10%. The distinction between BHAR as measured by 
Loin, and black hole growth rate as measured by Lom/LEdd can 
be confusing. Throughout the paper we refer to these quantities as 
"accretion rate" and "growth rate" respectively. 

In order to measure "growth rates" we estimate the black 
hole masses from the stellar velocity dispersion of the galaxies 
(a) by assuming the empirically calibrated Mbh — cr relation 
( Tremaine et al.ll2002h . We note that elliptical galaxies, classical 
bulges (bulges in galaxies with disks) and pseudo-bulges follow 



2 LEdd/L© = 3.3 x 1O 4 Mbh/M0 is the luminosity at which radiation 
pressure due to electron scattering balances the inward pull of gravity for 
optically thin spherical accretion onto the black hole. 



slightly different a — Mbuige relations, and therefore Mbh esti- 
mated from a Mbh — cr relation will b e diffe r ent from that estimated 
from a Mbh — Mbui ge re lation (Hu 2008; Gadotti & Kauffmann 
120091 : iGraham & Lill2009h . The offset can be attributed primarily 
to the presence of bars. However, as our results are based upon an 
average over a large sample of AGN, we do not believe this will 
introduce large biases. It is possible that bars are more prevalent in 
the youngest starbursts, in which case this effect will lead to slightly 
enhanced black hole mass estimates. We will note the implications 
of this effect where necessary in the results and discussions sec- 
tions. 

We separate our starburst galaxies into those wi th and without 
obscured AGN using optical emission line ratios (IB aid win et al.1 
1 19811 : iKauffmann etafl 120031 : iKewlev et all l2006h . For galaxies 
with a contribution to their emission lines from both star formation 
and AGN ("composite AGN" - the majority of AGN in our sample) 
we have removed the contribution to [O III] from star formation 
using an empirically calibrated mixing algorithm. In the following 
subsections we give further details of these procedures. Firstly, we 
justify the use of [O III] as an indicator of black hole accretion rate. 



3.1 [O in] as a black hole accretion rate indicator 

For obscured (Type 2) AGN observed in the optical wavelength 
range the high ionisation [O III]A5007 emission line is the best 
available measure of AGN power. The luminosity of this line, 
which originates in the Narrow Line Region, has been shown to cor- 
relate with numerous indicators of total AGN power in unobscured 
(Type 1) samples, where the central en gine can be viewed dir ectly 
to obtain bolometric luminosities (e.g. Mulcha ev et al For 
obscured (Type 2) AGN, [O III] luminosity correlates well with 
hard X-ray luminosity for Compton thin AGN, with a scatter of 
0.51d ex over 4 orders of magnitude in luminosity (Hec kman et al.1 
2005). Furthermore, because the AGN Narrow Line Region lies 
100's of parsecs from the central engine, [O III] is not affected by 
dust obscuration from the torus, unlike X-ray emission which can 
be heavily absorbed (Compton thick AGN). 

The conversion of [O III] luminosity into a bolometric AGN 
luminosity, from which we estimate black hole accretion rates, 
is based on the strong linear correlation between narrow [O III] 
emission and optical c ontinuum luminosity in unobscured AGN 
(Zakamska et al. |l2003l) . The conversion is discussed in detail in 
Hec kman et all (120041) and IKauffmann & Heckmanl (120091) . How- 
ever, we note here the implicit assumptions in this conversion for 
completeness. Firstly, the conversion is necessarily calibrated on 
unobscured AGN. By assuming this calibration holds for obscured 
AGN, we assume the AGN standard model in which [O III] is 
isotropically emitted and the difference between obscured and un- 
obscured AGN is in the orientation of the dusty torus relative to our 
sightline. Secondly, SDSS AGN samples extend to lower luminosi- 
tie s than the samp l es use d for the calibration. However, as shown 
by Heckman et al and discussed below, the total [O III] lu- 

minosity of a volume limited sample of AGN is dominated by AGN 
with high [O III] luminosities. Removing galaxies with [O III] lu- 
minosities below the limit to which the conversion has been verified 
has little effect on volume averaged total black hole accretion rates. 
Finally, it should be noted that at some level Loin should depend 
on the availability of gas clouds to ionise, and may therefore cor- 
relate with host galaxy properties. This may plausibly account for 
some of the scatter in the observed relations between Loin and 
other indicators of AGN strength. We do not believe it is likely that 
there is a substantial increase in the number of narrow line clouds 
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over timescales of several hundred Myrs, which could mimic the 
results found in this paper, especially given the relatively slow de- 
cline in star formation which presumably originates from the same 
gas supply. 



3.2 Missing lower luminosity and unobscured AGN 

There are two main ways in which the SDSS narrow line AGN sam- 
ple is an incomplete census of the local AGN population: (1) loss 
of low luminosity/low growth rate AGN when the lines become in- 
distinguishable from nebular emission from H II regions, or against 
the background stellar light; (2) broad line (unobscured) AGN are 
not included, because the bright nuclear light would preclude stud- 
ies of the integrated galaxy stellar light. 

Obscured, low luminosit y AGN con stitute the majority of the 
low-z AGN population jHao et al.ll2005l) . and recent deep X-ray 
surveys of local galaxies ar e revealing increa singly high fractions 
of low accretion rate AGN dGallo et al.ll2008b . Hard X-ray obser- 
vations are the most sensitive probe for (Compton thin) obscured 
AGN, and such low accretion rate systems would not be detected 
against the background stellar light in SDSS spectra. But, although 
these ultra low luminosity AGN are numerous, as noted above their 
contribution to the vol ume averaged growth of black holes is small 
(iHeckman et al.ll2004) . 

The volume averaged growth of black holes is do minated by 
high growth rate (L/LEdd) systems (see discussion in lHao et al.1 
2005). More specifically, Heckman et al. (2004) find that the low-z 
volume averaged black hole growth rate peaks for black holes with 
masses of ~ 10 7 ' 5 M© (their Fig. 1) and is dominated by systems 
with L/LEdd^ a few percent (their Fig. 3). Therefore, while miss- 
ing low-luminosity obscured AGN will have a negligible effect on 
estimated global BHARs, we may legitimately be concerned about 
missing high-luminosity unobscured AGN. Although the standard 
model for AGN suggests we can simply scale our results by the 
observed ratio of unobscured to obscured AGN, deviations from 
this model may affect our conclusions. We will return to this point 
where necessary in the discussion. 



3.3 Emission line analysis 

Emission line fluxes and errors, corrected for stellar contin- 
uum absorption, are extracted from the SDSS-MPA catalog 
(http://www.mpa-garching.mpg.de/SDSS/), see Brinchman et al. 
(2004) and Tremonti et al. (2004) for details of their measure- 
ment. All lines are corrected for dust attenuation using the observed 
Ha to H/3 ratio and the two component dust attenuation curve of 
Chariot & Fall (2000), as presented in Wild et al. (2007) and de 
Cunha et al. (2008). For star forming galaxies and composite- AGN 
(where emission from both star formation and an AGN is present), 
an intrinsic Ha/H/3 ratio of 2.87 is assumed; for pure- AGN we use 
a slig htly higher Ha/H/3 = 3.1 (IOsterbrockl[T98^ : iKewlev et al.l 
2006). A very small fraction of our starburst galaxies have lines 
that are too weak to be measured accurately, about 4% of the star- 
burst sample with ages younger than 600 Myr and 0.5% younger 
than 250 Myr. These fractions are too low to cause any bias in our 
results. 

In Fig.[3]we present the [ N HlA6585/Ha vs. [O m]/H/3 line ra- 
tio diagram (the BPT diagram, Bald win etalTl 98D for all galaxies 
in our sample with lines measured at > 3a confidence. This dia- 
gram is commo nly used to separate purely star formi ng galaxies 
from AGN (e.g. IKewlev et allboOll : iKauffmann etai1 l2003). The 




Figure 3. The BPT diagram. Bulge-galaxies with all four emission lines 
measured at > 3cr are shown as gray-scale. The coloured circles are the 
starburst sample, with colour indicating starburst age as in Fig. 1. The lines 
indicate the positions of the demarcation lines Csf (dashed red), Cagn 
(dashed blue), CsF/comp (lower black) and C CO mp/AGiV (upper black) 
as described in the text. The white dotted lines show how we assign k to 
composite- AGN when the AGN and starforming contributions to the emis- 
sion lines are calculated (Eqn.0. 
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Figure 4. An important consequence of our results for the study of the 
AGN- starburst connection in individual objects using wide aperture inte- 
grated optical spectra alone, is that "pure" AGN will not be found in star- 
burst galaxies, but only appear after the starburst. Before this time the star- 
formation will cause AGN to manifest as composite galaxies. Here we show 
the number of starforming (blue), composite- AGN (orange) and pure- AGN 
(red) in our sample as a function of time since the onset of the starburst. 



starburst sample is overplotted, with colour indicating starburst age. 
We note that the trend for the youngest starbursts to have lower 
metallicities evident in this Figure is expected based on previous 
studies and simulations (Kewley, Geller & Barton 2006; Rupke, 
Kewley & Barnes 2010). 

The upper black line lies close to the theoretical "maximum 
starburst" contribution from Kewley et al. (2001), above which 
Seyferts (high [O lll]/H/3) and LINERS (low [O lll]/H/3) are eas- 
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Figure 5. For composite- AGN we must separate the contributions to [O in] 
and Ha emission from AGN and star formation. On the left is the distribu- 
tion of total [O in] vs. Ho; luminosities for AGN (red), star forming galaxies 
(blue) and composite- AGN (orange). From outside to inside the contours 
encompass 90, 60, 30 and 10% of the data points. On the right, [O in] and 
Ho; luminosities for the composite- AGN have had the contribution from star 
formation removed as described in the text. We see that this correction pro- 
cedure places the composite- AGN in the same range of [O in] luminosities 
as the pure AGN. 



ily identified as plumes. As shown by Kauffmann et al. (2003) the 
ordinary star forming galaxies in the SDSS define a tight sequence 
that lies far below the theoretical model maximum, close to the 
red-dashed line indicated in the Figure. The simplest interpreta- 
tion of the objects which lie between these two populations is that 
both starformation and AGN contribute to their emission lines, so- 
called "composite AGN". This is supported by their younger ages 
as measured from the stellar continuum. The black lines indicate 
our separation between pure star forming galaxies and composite 
AGN (lower black line), and composite AGN and pure AGN (upper 
black line). The blue-dashed line is placed to cross through the cen- 
ter of the pure- AGN. The red-dashed line is placed to cross through 
the center of the star-forming galaxies 

In this work, all demarcation lines are defined using the same 
equation for simplicity: 

(i) 



BPT, 



1.3 + 



BPTcc - C 

where the four lines in Fig. [3] are defined by Csf = —0.08 
(red), CsF/comp = -0.03 (lower black), C comp/A GN = 0.5 
(upper black), and Cagn = 0.9 (blue). We place the star- 
forming/composite demarcation line a little lower than that used by 
Kauffmann et al. for their sample. We will discuss the implications 
of the placement of this line on our results in the next subsection. 

Galaxies which lie below Csf /com P have negligible contribu- 
tion to their [O III] line from an AGN and we classify them as star 
forming galaxies, and those that lie above C comp /AGN have neg- 
ligible contribution from star formation and we classify them as 
pure- AGN. Those that lie in between these two lines we classify 
as composite AGN. In Fig. [4] we plot the number of each emission 
line class vs. the age of the starburst, from which it is clear that in 
optical spectra probing a large spatial area pure- AGN will only be 
observed about 500 Myr after the onset of the starburst due to con- 
tamination of the AGN lines by the strong nebula emission from 
star formation. Before this time, AGN will necessarily be classified 
as "composite- AGN". 

In this paper, we include both composite and pure AGN in 
our measurements of the total BHAR and distributions of L/LEdd- 
While galaxies which lie below the CsF/com P demarcation line 
may contain AGN, these will have very low L/LEdd, therefore this 
line should be thought of as a limit on L/LEdd of our AGN sam- 
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Figure 6. In the left panels we show the AGN contribution to Loin 
recovered from our Monte Carlo simulations combining AGN and star- 
formation-dominated galaxies using the correction method described in the 
text, compared to the true AGN Loin, for two different ranges of starburst 
age. In the right panels, we show the position of all the simulated composite- 
AGN on the BPT diagram, with the same demarcation lines as shown in Fig. 
[3] The points are colour coded by the AGN contribution to the total Lqiii 
of the simulation galaxy. 



pie, rather than a firm division between galaxies with and with- 
out an accreting black hole ( see also the detailed discussion in 
Kauffmann & Heckman 2009). We will return to this point later 
when we account for the incompleteness of our AGN sample as a 
function of L/Lsdd- 

3. 3. 1 Correction ofLoiu for contamination from starformation 
in composite AGN 

In order to measure the AGN luminosity of composite AGN from 
the [O III] emission line, we must correct for contamination due to 
starformation. Likewise, to measure star formation rates from Ha 
luminosities, we must correct for contamination from the AGN. 
Our method is based upon the assumption that as AGN contribu- 
tion increases, galaxies move diagonally away from the starform- 
ing sequence approaching the AGN demarcation line (following the 
white dotted lines in Fig. [3}. 

Inverting Eqn. 1 to obtain C for an individual composite AGN, 
we measure its distance from the starforming sequence to be: 

C — Csf 



Dbpt 



(2) 



Cagn — Csf 

which is normalised such that galaxies that lie in the center of the 
star-formation- dominated bulges (red dashed line) and in the center 
of the AGN-dominated bulges (blue dashed line) have Dbpt = 
and .Dbpt = 1 respectively. 

Observing that both pure AGN and star forming galaxies in 
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our sample show very tight relations between Ha and [O III] lumi- 
nosities (Fig. [5] left panel), and making the assumption that Dbpt 
defined in this manner is proportional to the fraction of Ha lumi- 
nosity arising from the AGN, we derive empirically the following 
correction for [O III] and Ha in the composite AGN: 

£oiII,AGN = £oiII,Tot — — -DbPT )Z^Ha;,Tot (3) 

^Hq,SF — ^Ha,Tot — ^BPT^Ha,Tot (4) 

where Tot indicates the total contribution to the emission luminos- 
ity from both AGN and star formation, k is the [O Ill]/Ha ratio 
of starforming galaxies at the point on the starforming sequence 
where the composite AGN originated from (following the dashed 
white lines to meet the red dashed line in Fig. [3]) 



log 10 (2.87^) 



BPT* 



= 1.3 + 



0.61 



BPT X SF - Csf 



BPTy + m(BPT x ~ BPT X ) 



(5) 
(6) 
(7) 



which we solve numerically, setting m — 1.3 and as suming an 
intrinsic Ha to H/3 flux ratio of 2.87 ( Osterbrock 1989). The suc- 
cess of the separation procedure can be seen in the right panel of 
Fig. \5\ where the corrected [O III] to Ha line luminosities of the 
composite- AGN follow the relation of the AGN. 



3.3.2 Monte Carlo simulations of the Loin correction procedure 

Because the correction of Loin for contamination from star forma- 
tion is not yet common practice in the literature, and has not been 
thoroughly tested, we present here a simulation to demonstrate the 
validity of the technique. We proceed by taking all pure- AGN with 
starburst ages > 250Myr, and add their lines one-by-one into each 
pure star forming galaxy in our starburst sample. We do this for dif- 
ferent age bins of the star-formation-dominated starburst galaxies 
and repeat the experiment for a control sample of blue- sequence 
bulge-galaxies. This control sample is extracted from the region 
defined by the dotted box in Fig. [T] and is matched to the starburst 
sample in black hole mass, stellar mass, redshift and stellar surface 
mass density. We then recalculate the position of the galaxies on the 
BPT diagram, perform our [O III] separation analysis and compare 
the output AGN component of the total [O III] luminosity, to the 
input true AGN luminosity. 

The result is shown in Fig. [6| for two starburst samples with 
age bins of < age/Myr < 200 and 400 < age/Myr < 600, 
and the control sample of blue- sequence bulge-galaxies. Similar 
results are obtained for the intermediate age starburst sample. The 
left hand panels show that on average we recover the input Loin 
well. The right hand panels show the simulated composite-AGN 
on the BPT diagram, colour coded by the fraction of Loin con- 
tributed by the AGN. On average we see that, as expected, galaxies 
close to the SF-composite demarcation line (lower black line) have 
small contributions, whereas those that sit close to the composite- 
AGN demarcation line (upper black line) have contributions close 
to unity. 

It is noticeable that the youngest starbursts suffer from an 
overestimation of Loiii,agn by about a factor of 2 on average. 
This results from their lower mean metallicity which, if mixed with 
a LINER ratio AGN, means they do not follow the assumed mix- 
ing tracks on the BPT diagram. When mixed with Seyfert line 
ratios, however, the method works accurately. We emphasise that 
the method presented here is suitable for the majority of bulge- 
dominated galaxies at low-redshift. However, for low-metallicity 




0<Age/Myr<250 
250<Age/Myr<600 
Blue sequence 



0.010 0.100 
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Figure 7. Result of the Monte Carlo simulation to estimate the fraction of 
AGN recovered as a function of normalised BHAR (L/Lgdd) an d starburst 
age. The dotted and dashed lines are for the starburst sample with ages as 
indicated in the figure. The solid line shows the result for the control sample 
of ordinary bulge dominated starforming galaxies. By 250 Myr after the 
starburst, the fraction of AGN recovered is similar to the control sample, 
however, at younger ages we can recover fewer AGN due to the masking of 
AGN lines by the emission lines from starforming regions. In our analysis 
we only include bins which are >50% complete. 



star-forming regions there is ultimately a degeneracy that prevents a 
unique solution for the fraction of Loin originating from the AGN 
using these line ratios alone. In Section |4] we deliberately exclude 
the very youngest starbursts (<30Myr) from some parts of our anal- 
ysis, which reduces the impact of this degeneracy as these galaxies 
have the lowest metallicities (see Fig. [3] and Section [33l >. We note 
in the results and discussion sections where this may affect our re- 
sults. 



3.4 AGN Completeness 

There are three possible causes for not identifying AGN in a purely 
optical sample: (1) the AGN emission lines are completely con- 
cealed by the strong nebular lines from star formation; (2) the AGN 
emission lines are not identifiable against strong stellar continuum 
radiation; (3) the AGN emission lines are obscured by dust. For the 
results presented in this work, we are primarily concerned about 
correcting for the first effect, which causes us to preferentially lose 
weaker AGN in the youngest starbursts. Effect (2) will not cause a 
similar bias as the red broad-band magnitudes of the sample change 
little with the time probed and the low AGN luminosities involved 
do not contribute significantly to the time averaged BHAR as dis- 
cussed in Section [3Jl Effect (3) may cause us to loose a small num- 
ber of objects at older ages (> 350Myr) corresponding to the onset 
of the Asymptotic Giant Branch phase in lower mass stars, where 
we start to see a significant increase in the mean B aimer decrement 
of the starburst galaxies. This trend of dust content with starburst 
age is the subject of future work. 

Returning to the first effect, we use the simulations described 
above to calculate the completeness of our samples as a function 
of L/LEdd by calculating the position of the simulated composite 
AGN on the BPT diagram and finding whether they would be in- 
cluded in our AGN sample or not (i.e. lie above the lower black line 
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Figure 8. Ha luminosity from star formation, corrected for dust attenuation 
and contamination from AGN emission as described in the text. Blue points 
indicate star-forming galaxies, orange indicate starforming-AGN compos- 
ites. Upper limits indicate unclassified objects (black) and pure- AGN (red). 
The gray solid line and dotted gray lines trace the 50th, 10th and 90th per- 
centiles of the distribution using a running bin width of 71 galaxies. The 
black dashed line shows an exponentially declining star formation rate with 
r = 0.3 Gyr, normalised to the median SFR in the first 40Myr. Due to 
the significant fraction of objects at late ages for which we only have upper 
limits of L(Hq.) the 10th percentile is only indicated to 400 Myr. 



of order 10-30 Myr, an almost flat "coasting" phase which lasts 
to around 400 Myr, and a subsequent decline. Throughout this pa- 
per, we have parameterised the overall decline in starformation by 
an exponential. In this Figure an exponential with a timescale of 
t — 0.3 Gyr is overplotted as a black dashed line. It is normalised 
to the median Ha luminosity in the first 40 Myr. As discussed in 
Section when fitting starburst ages to our bulge-galaxies we 
assume a range of r values for the starburst models. This range 
was decided upon by taking the observed decline in Ha luminosi- 
ties into account: allowing r values much shorter (longer) causes 
the observed Ha luminosities of the older starbursts to lie above 
(below) the predictions of the models used to estimate the starburst 
age. While neither stellar continuum nor line luminosities alone 
constrain the fall in star formation during a starburst, it is encour- 
aging that a simple model shows that they can be consistent with 
each other. A more precise study of the decline in star formation 
following a starburst will be feasible with the data and methods 
presented here. 

The starburst timescales of several hundred Myrs measured 
here are consistent with measurements from two entirely indepen- 
dent m ethods. Firstly, from obse r ving s ta rbursts in close galaxy 
pairs (iBarton. Geller. & Kenvonl l2000l : iFreedman Woods et all 
2010), and secondly fro m the stellar mass s urface density profiles 
of elliptical galaxies (iHopkins & Hernquisll2009h . In the latter pa- 
per a very similar form for the star formation history is also inferred 
for the component of elliptical galaxies thought to be formed during 
gas rich starbursts. 



in Fig. O. The resulting measured completenesses, shown in Fig. 
13 are then included in all quoted values of the AGN fraction. We 
can see that even at the youngest starburst ages we recover >90% 
of AGN with L/L E dd > 0.01, and at later times this rises to 95%. 
Therefore, we are not losing a significant number of high accretion 
rate AGN from our samples. However, it is clear that AGN with 
L/LEdd < 0.003 are difficult to detect at all ages during the star- 
burst, with recovery rates of <50%. In order to investigate further 
the appearance of low accretion rate AGN, follow-up observations 
at higher spatial resolution and, possibly, in other wavebands will 
be necessary. 

In order to be sure that our key results are completely robust to 
changes in the precise positioning of the demarcation between star 
forming galaxies and AGN (lower black line in Fig. [3]), we repeat 
the key results but this time allowing all galaxies which lie above 
the red line in Fig.[3]to be composite- AGN. We find that although, 
of course, there are many more "AGN" in our sample these have 
low growth rates and the total BHAR remains the same. 

3.5 Decline in SFR of a starburst 

Not only does our starburst sample allow us to investigate the time 
averaged AGN activity during a starburst, as presented in the re- 
sults section below, it also provides a means by which to measure 
the form of the decay in star formation after the onset of a star- 
burst. Very few constraints currently exist on this decay timescale, 
and instantaneous or very short starbursts are often assumed in the 
literature which are not consistent with the results presented here. 

In Fig. [8] we present the Ha luminosity of all starforming and 
composite galaxies in our starburst sample, as a function of star- 
burst age derived from the stellar continuum. Composite galaxies 
have had the contribution to Ha luminosity from AGN removed 
(Eqn.|4]). We can clearly identify 3 phases: a short "starburst" phase 



4 RESULTS 

By averaging over a sample of 400 bulge-galaxies, we can build a 
picture of the time-averaged BHAR with respect to the global star 
formation activity of the sample, even though individual accretion 
events may be relatively short lived. To calculate the average black 
hole accretion rate occurring in each time bin, we sum the total 
Loin originating from the pure- and composite- AGN corrected for 
contamination from star formation, and divide by the total number 
of starburst galaxies in that bin, i.e. including the pure starbursts 
currently not experiencing a strong black hole accretion event. In 
the top panel of Fig. [9] we plot the mean BHAR (M /yr/galaxy) as 
a function of time since the onset of the starburst, finding it to re- 
main low during the early phase of the starburst, and increase only 
at ages >250 Myr. A similar offset between starburst and AGN ac- 
tivity has been sugg ested for a handful of nearby galaxies studied 
in detail (e.g iTadhunter et al.ll996l : lEmonts et al.l200rj|Gotd2006l ; 
iDavies et~al] 20071) . In the lower panel of Fig. [9] we show how the 
ratio of SFR/BHAR evolves rapidly with time since the starburst, 
where SFR is measured directly from Ha luminosity. We find that 
the total mass of stars formed is 2 x 10 3 times the total mass ac- 
creted onto the black hole during the first 600 Myr after the star- 
burst. This ratio is a little higher than, but of the same order as 
the ratio of stellar mass to black hole mass in present-day bulges 
(Haring & Rix 2004), suggesting that the processes we observe in 
the starbursting bulges may be more broadly relevant to the co- 
evolution of bulges and black holes. We will return to this point in 
the discussion. We have deliberately not plotted the average BHAR 
for starbursts with ages younger than 30 Myr in Fig. [9] (20 galax- 
ies), because of the additional incompleteness for these galaxies 
with extreme star formation rates and the additional uncertainty in 
the separation of [O III] into starformation and AGN components 
(i.e. the initial spike of galaxies with strong Ha fluxes in Fig.[8}. 
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Figure 9. Top Panel - the mean BHAR averaged over all galaxies in the sam- 
ple (black line), estimated from a robust mean with 5a clipping. The gray 
shaded area shows the typical 10th and 90th percentile range on the total 
BHAR estimated from bootstrap resampling of the data. The initial 30 Myr 
after the onset of the starburst is not plotted, as explained in the text. The 
model prediction for the average gas mass loss rate from the stars formed 
during the starburst (solid gray line) assumes an accretion efficiency onto 
the black hole of 0.5% of the ejecta from low-mass stars between the age 
of 250 and 600 Myr. Stars return mass to the interstellar medium through 
SNe explosions and fast winds from massive stars (dotted line) and plane- 
tary nebula ejections and stellar winds from lower mass stars (dashed line). 
Bottom panel - the ratio of SFR (from Ha) to BHAR, with mean value and 
errors calculated as for the top panel. 
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Figure 10. The AGN fraction in bulge-galaxies with starburst ages 
<250Myr (dotted line), with starburst ages between 250 and 600 Myr 
(dashed line) and for bulges undergoing ordinary star formation and 
matched to the starburst sample in stellar mass, surface mass density, black 
hole mass and redshift (solid line, see dotted box in Fig. [I}. Left: As a func- 
tion of BHAR. Right: As a function of black hole growth rate (L/LEdd)- 
In both panels the AGN fractions have been corrected for incompleteness 
caused by the concealment of AGN emission by emission from star forming 
regions, and only bins which are >50% complete are plotted. 



|7] and Section 13.41 ) and only bins with an overall completeness of 
greater than 50% are plotted. Comparing the distribution of BHAR 
(left panel) in the three samples we see that both young and old star- 
bursts show an enhanced number of powerful AGN with respect to 
the control blue- sequence sample, which has been matched to the 
starburst sample in redshift, stellar mass, stellar mass surface den- 
sity and black hole mass. However, in the youngest starbursts the 
AGN fraction is reduced at all accretion rates relative to the AGN 
fraction in the older starbursts. In both young and old starburst sam- 
ples the distribution of BHAR turns over at low BHAR, suggesting 
that low accretion rate AGN in starbursts are suppressed compared 
to in ordinary star forming bulges. However, given the increasing 
importance of completeness corrections at lower BHAR this result 
should be verified with additional observations. 

The right hand panel of Fig. [10] provides a different view of 
black hole accretion patterns during the starburst. Here we see 
that the fraction of starbursts with rapidly growing black holes is 
enhanced relative to our control sample of blue- sequence bulge- 
galaxies. Specifically, the fraction of galaxies with black holes ac- 
creting at rates greater than 0.7% of the Eddington rate is ~ 1.2 
times higher in starbursts with ages between and 250 Myr, than 
in ordinary blue-sequence bulge-galaxies, and ~ 2 times higher in 
starbursts with ages between 250 and 600 Myr. Although incom- 
pleteness prevents us from probing to very low growth rates, the 
turn-over in the distribution of L/LEdd occurs at higher L/Lsdd in 
the youngest starbursts, compared to both the older starbursts and 
the control sample. The shape of these distributions strongly sug- 
gest that low growth-rate accretion events are suppressed during the 
starburst. 

We have not included the 7 composite-AGN with ages 
younger than 30 Myr in the dotted histograms of Fig. [Toj for the 
reasons explained above. The effect of including them is small, and 
does not change the impression given by the figures. As discussed 
in Section[3] it is possible that the presence of bars at certain phases 
during the starburst may cause biases in the right panel of Fig. [Toj 
due to our estimation of black hole mass from stellar velocity dis- 
persion. However, the most likely consequence is that we slightly 
overestimate the black hole mass of a small fraction of galaxies 
with the youngest starburst ages. Correcting for this would only 
enhance our conclusion that the highest growth rate events are not 
suppressed in young starbursts. 

Combining all these results we conclude that the global ac- 
cretion onto black holes is lower in young starbursts than in older 
starbursts. This is caused by an overall drop by around a factor of 
two in the number of accretion events at all BHARs above at least 
10 _3 M©/yr. The possible trend for the turnover in the distribution 
of L/LEdd to occur at higher values of L/LEdd in starbursts com- 
pared to the control sample, and at higher values in younger com- 
pared to older starbursts, should help constrain the physical mech- 
anism responsible for this decrease. We will return to this point in 
the discussion. 



To uncover the primary factor driving the low global accretion 
at young starburst ages, in Fig.[l0]we plot the distribution of BHAR 
(left) and growth rate (right) in the starburst and control samples. 
We have divided the starburst sample into two at 250 Myr and we 
compare to galaxies with more ordinary recent star formation his- 
tories (blue- sequence bulge-galaxies matched in black hole mass, 
stellar mass, redshift and stellar surface mass density, see dotted 
box in Fig.[T]). All histograms have been corrected for incomplete- 
ness in the samples as a function of L/LEdd and starburst age (Fig. 



4.1 Summary of observational results 

The existence of some connection between starbursts and AGN has 
been debated for several decades. Sanders et al. (1988) argue for an 
evolutionary connection on galaxy wide scales between ULIRGs 
and QSOs. Heckman et al. (1997) argue for a causal connection on 
nuclear scales (Gonzalez Delgado et al. 2001; see Veilleux 2001 for 
a review of the subject). Here we argue for its existence for the first 
time based on a large statistical sample with a well defined control 
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sample. However, we modify the nature of connection in several 
ways. To summarise the key observational results, we find that: 

• Galaxies with bulges experiencing a strong starburst are more 
likely to host a rapidly growing and rapidly accreting black hole 
than galaxies with bulges experiencing more normal levels of star 
formation. 

• The total mass of stars formed in the first 600 Myr after the 
onset of the starburst is 2 x 10 3 times the total mass accreted onto 
the black hole. The black holes increase their mass by on-average 
5% during this time. 

• On average, the rate of accretion onto black holes rises steeply 
about 250 Myr after the onset of a burst of star formation. 

• The lower average BHAR in young starbursts compared to old 
starbursts is caused by a decrease by about a factor of two in the 
number of accretion events at all BHARs above at least lO 3 M /yr. 

• Alternatively, our results can be viewed as an enhancement in 
the global accretion onto black holes during a starburst, relative to 
that occurring in ordinary blue-sequence bulge-galaxies, but mainly 
in bulges with starbursts older than 250 Myr. 

• On the rare occasions that a black hole accretes during the 
early phase of the starburst, it grows rapidly (has L/LEdd> 0.01). 

• While higher growth rate events are enhanced compared to 
blue-sequence bulge-galaxies, lower growth rate events instead ap- 
pear suppressed. 

We note that our results apply primarily only to black holes in the 
mass range 10 6 5 — 10 7 ' 5 M©, moderate luminosity AGN (bolomet- 
ric luminosities of order 10 42 — 10 44 erg/s), and powerful starbursts 
which increase the stellar mass in the central <2kpc radius of the 
galaxy by at least 10%. 



5 DISCUSSION 

5.1 Unobscured AGN 

To begin our discussion, we return to one of the assumptions im- 
plicit in our calculations: to extrapolate our results for obscured 
(Type 2) AGN to the growth of all black holes, the simplest ap- 
proach is to assume the standard unified model for AGN in which 
the difference between unobscured and obscured AGN is in the 
viewing angle alone. Under this assumption, the effect of our ex- 
clusion of unobscured (Type 1) AGN will mean that we have un- 
dercounted the total amount of black hole growth in our population 
of starburst galaxies. 

There have been many estimates of the ratio of Type 2 to Type 
1 AGN in the local universe, with values ranging between ^3 (e.g. 
IdeGriip etallll992h and -1 (e.g. Hao et al. 2005). Thus, we would 
need to increase the rate of black hole accretion by factors rang- 
ing from 1.3 to 2. There is some dispute about whether the relative 
number of Type 1 and Type 2 AGN is an increasing function of 
luminosity at low redshift. For example, de Grijp & Miley (1992) 
found no evidence for such a dependence based on the mid-IR lu- 
minosity functions, while Hao et al. (2005) found that the ratio did 
increase by a factor of 2-4 over a range of about four orders of mag- 
nitude in [O III] luminosity. Thus, we may be underestimating the 
relative importance of high accretion rate events by examining only 
obscured AGN. 

A more interesting possibility would be that the Type 1 and 
Type 2 AGN are related to one another in an evolutionary sense 
(and not distinguished purely by viewing angle). For example, it 
is physically plausible that obscured AGN evolve into unobscured 



AGN as feedback associated with the AGN and/or the surrounding 
starburst clear away much of the obscuring material (e.g. Sanders 
et al. 1988). In this case, by excluding Type 1 AGN we may be 
selectively undercounting the relative amount of black hole growth 
during the late phases of the starburst compared to the early stages. 
If so, and if this scenario were applicable to the moderate luminos- 
ity AGN studied in this paper, then this would only increase the 
sense of the increase in black hole growt h at these late times . 

It is interesting in this regard that iDavies et ID (2007) have 
measured the nuclear starburst ages of a sample of mostly unob- 
scured local AGN, finding a time delay between the onset of the 
starburst and the onset of significant black hole accretion. The mag- 
nitude of the delay was similar, although a little short e r than , that 
found in our study. As pointed out by IDavies et aD ([2007| ), the 
Galactic center also fits into this picture, where the winds from 
nearby young OB stars are not being accreted as efficiently as ex- 
pected (Baganoff et aD[2003; Paumard et al. 20q3). 



5.2 Black hole accretion from stellar mass loss 

One popular hypothesis for the physical mechanism responsible for 
the feeding of black holes is through th e accretion of matter ejecte d 
during stellar evolutionary processes dNorman & Scovilldll988r) . 
We can investigate the link between starburst activity and BHAR by 
examining the mass returned to the interstellar medium (ISM) by 
supernova (SNe) explosions, planetary nebula ejections and stellar 
winds predicted from stan dard stellar population synthesis models 
(iBruzual & Charlotll2003l Chariot & Bruzual, in prep.). The solid 
gray line in Fig. [9] shows the mass loss rate of stars formed during 
a starburst, and the dotted and dashed lines separate this into mass 
loss from high mass (predominantly SNe) and lower mass starjf|. 
During the initial phase of the starburst mass loss from SNe ex- 
plosions and fast winds from O/B stars dominates the mass loss 
budget, however as can be seen in Fig. [9] mass loss from lower 
mass stars soon becomes dominant. The model is normalised to 
match the total mass of gas per galaxy converted into stars in our 
starburst sample, which we have measured f rom the Ha lin e lu- 
minosity, converted into star formation rate feennicuttlll998h and 
corrected for contamination from AGN. The star formation rate of 
the model declines exponentially with a timescale of 300 Myr, in 
agreement with the observed fall-off in the Ha luminosity arising 
from star formation in the sample. 

The model for stellar mass loss rate has been placed onto 
the same scale as the BHAR by assuming that 0.5% of the mass 
returned to the ISM by evolved low and intermediate mass stars 
within 250 and 600 Myr is accreted by the black hole. This fraction 
is similar to that predicted by, for example, the model of Ciotti & 
O striker (2007) and to that inf erred in galaxy bulges with pre dom- 
inantly old stellar populations ( Kauffmann & Heckman 2009). It is 
immediately evident that the average BHAR of the sample does not 
follow the total mass loss rate from stellar processes (solid line), 
in particular during the early period of rapid mass loss from SNe 
(dotted line). One key difference between mass loss from SNe com- 
pared to that from lower mass stars is in the ejecta velocity: several 
thousand km/s for SNe compared to several tens of km/s for lower 
mass stars. Fig.[9]suggests that the fast ejecta from SNe must escape 
the nuclear region without interacting with the central black hole. 



3 Mass loss from low mass stars in the pre-existing old stellar population 
in the bulge is not significant during these phases of a strong starburst. 
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This is consistent with the model of Norman & Scoville (1988) and 
with detailed observations of nearby AGN (Davies et al. 2007). 

If a simple model for the feeding of black holes through mass 
loss from slow stellar ejecta is correct, then an important question 
is whether SNe are simply inconsequential for black hole growth, 
or whether they cause enough disruption of the ISM to prevent the 
accretion of the slower stellar ejecta by the black hole. Two features 
of our results suggest that the latter may indeed be the case. Firstly, 
from Fig. [9] we see that the BHAR does not track the mass loss from 
intermediate and low mass stars at the earliest times between 100 
and 200 Myr after the starburst. During this time, supernova rates 
remain high and may provide the energy required to inhibit the ac- 
cretion of slower stellar ejecta and any residual gas. Secondly, in 
the right panel of Fig.[l0]we see that the distribution of AGN frac- 
tion as a function of black hole growth rate changes with starburst 
age, i.e. while higher growth rate events are enhanced compared 
to blue-sequence bulge-galaxies, lower growth rate events are in- 
stead suppressed. This pattern is stronger in young starbursts than 
in old starbursts, and the turn-over in the distribution of L/LEdd oc- 
curs at higher values for young starbursts than for old starbursts. A 
scenario in which growth rate is linked to the mass of an accreted 
cloud, and supernovae preferentially disrupt/expell the lowest mass 
clouds, would account for the observed suppression of low growth 
rate events during starbursts. 

5.3 Alternative models for black hole accretion 

Our observations have direct implications for theoretical models 
that attempt to explain t he joint evolution of black holes and galax- 
ies dSilk & Reesl l l998l; lHaehnelt et all 1 19981 : iHopkins et allbood : 
ICiotti & O striker 20q3). Our results clearly show that the rate of 
growth of the black hole is offset in time with respect to the total 
rate of stellar mass-loss. Moreover, the ratio of the mass accreted by 
the black hole to the mass lost by intermediate mass stars alone rises 
as the starburst ages. However, there are no direct observational 
constraints on the primary physical mechanism responsible for the 
feeding of a black hole. In the previous subsection we have chosen 
to focus our comparison on the scenario in which black holes feed 
primarily on stellar ejecta from intermediate-mass stars, due to its 
simple application to our data through stellar population synthesis 
modelling. However, we have found that this model in its simplest 
form is not sufficient to explain our results. A plausible additional 
component to the model, allowing it to match our results, is that su- 
pernovae disrupt the gas in the nuclear regions and prevent the early 
accretion of slower material ejected by lower mass stars. We now 
discuss two alternative scenarios which are not currently advanced 
enough to provide direct comparison with our data. 

One alternative for the origin of the time delay is that stellar 
ejecta are initially intercepted by gas in the nucleus not yet con- 
verted into stars. Numerical simulations which could test this sce- 
nario are not yet available. A second possibility is that the delay in 
black hole growth could be purely dynamical in nature, relating to 
the collision of two bulges during a merger. Some of the images in 
Figs. lAll and lA2l show galaxies which have experienced significant 
recent disruption. In simulations of major mergers with a black hole 
in each progenitor galaxy, rapid accretion occurs after coalescence 
of these black holes, which lag s behind the peak of the starburst by 
a sho rt time (~ O.lGyr, e.g. Matteo et al.ll2005l ; Uohansson et all 
2009). Although this delay appears too short to match our results, 
higher spatial resolution simulations should be developed before 
this possibility is ruled out. Simulations are needed that better re- 
solve the complex physical processes occurring in a galactic nu- 



cleus during a starburst. Although both of these alternative scenar- 
ios for the feeding of black holes may plausibly lead to a delay in 
accretion after a starburst, the intriguing suppression of low growth 
rate accretion events during the starburst may still require the pres- 
ence of some form of feedback process. In order to gather evidence 
for or against each of the scenarios discussed here, spatially re- 
solved spectroscopy would be invaluable, allowing us to identify 
double nuclei, measure starformation rates as a function of radius, 
and obtain more accurate accretion rates. 

Kauffmann & Heckman (2009) found that the distribution in 
the growth rates of black holes at low-redshift swaps from a power- 
law distribution for black holes in old bulges (little or no star- 
formation) to a log-normal distribution for black holes in young 
star-forming bulges. They showed that the accretion rates in the 
old systems were consistent with the black hole accreting 0.3 to 
1% of the mass loss from evolved stars in the bulge (consistent 
with the value we derive). They speculated that the transition to the 
log-normal regime was associated with the black hole regulating 
its own growth in the young gas-rich bulges. In this paper, we are 
suggesting that stars may not only provide the fuel source, but also 
the feedback (in the form of supernovae associated with actively 
star-forming bulges). 



5.4 Implications for Starbursts 

Our sample of starburst galaxies has other useful applications, be- 
sides tracing the growth of black holes, such as measuring the de- 
cay in star formation following a starburst (Section [33} . It is clear 
both from these observational results, and from simulations of star- 
bursts induced by mergers that the assumption of an exponential 
for the decay in star formation following a starburst is an oversim- 
plification. Given the apparent complexities in star formation histo- 
ries suggested by these simulations, it is perhaps surprising that our 
simple approach produces noticeable trends and easily interpretable 
results. With the development of galaxy simulations with higher 
spatial resolution, the inclusion of more realistic physical prescrip- 
tions of star formation, and the conversion of output into "observ- 
able" properties, in the near future we should be able to compare 
simulations directly with the sample presented in this paper to fur- 
ther constrain the true star formation histories of starbursts. 



5.5 Global implications 

Averaged over the full sample, we find that the black holes increase 
their mass by 5% during the first 600 Myr after the starburst. While 
clearly not the dominant phase of the growth of the black holes, 
this value is far from insignificant. We have purposefully restricted 
our sample to the strongest starbursts that occur in bulges in the 
local Universe, simply to enable an accurate measure of the star 
formation timescales. However, inspection of the images in Figs. 
IAll and lA2l implies that our sample is not composed exclusively of 
rare major mergers, and the same time delay may equally apply to 
more ordinary episodes of star formation. At higher redshift, where 
gas fractions are higher, and star formation levels are enhanced, the 
same patterns observed in the strongest starburst systems at low 
redshift are expected to become more widespread. This would im- 
ply an increased fraction of high accretion rate AGN with redshift 
and, possibly, a decrease in the fraction of low accretion rate AGN. 
Growth rates would similarly be enhanced. 

Over the lifetime of a galaxy the mass loss budget from stellar 
evolutionary processes is dominated by intermediate and low mass 
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stars, thus our results are consistent with a scenario in which the 
tight observed correlation between black hole and bulge mass is 
maintained through the consumption by the black hole of a fixed 
fraction of low-velocity winds expelled by such stars in the bulge. 

To evaluate this quantitatively, we adopt a simple prescription 
in which we follow the evolution of the stellar mass loss and black 
hole growth over a period of lOGyr (the age of typical present- 
day bulges) in a model exponentially declining star formation his- 
tory with decline timescale of 300 Myr. For simplicity, we assume 
that none of the supernovae ejecta are accreted by the black hole 
and that none (0.5%) of the low- velocity stellar ejecta are accreted 
before (after) a time of 250 Myr. The late-time accreted fraction 
we adopt i s consistent with observations of black hole growth in 
old bulges (Kauffmann & Heckman 2009). This simple model pre- 
dicts that at a time of lOGyr, the ratio of the remaining stellar 
mass to the accumulated black hole mass will be ~ 400. Obser- 
yations of present day bulges yield an actu al mean ratio of 700 
dMarconi & Hundl2003l : lHaring & Rixll2004l) . Given the uncertain- 
ties in converting our measured [O III] luminosities into black hole 
accretion rates, uncertainties in the stellar mass-loss rates, and the 
simple model we have adopted, the similarity of these two values is 
quite intriguing. It suggests that the processes at work in our local 
starburst sample may well be relevant to the co-evolution of black 
holes and bulges over cosmic time, and that black hole growth has 
been fueled predominantly by mass loss from intermediate mass 
stars. 

If stars are the primary source of fuel for growing the black 
hole and also play a major role in the feedback processes that limit 
this growth, then it starts to be less puzzling that black hole mass is 
linked to the stellar mass of the bulge. 
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APPENDIX A: SDSS IMAGES OF STARBURSTS 
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Figure Al. SDSS postage stamp images of the 20 youngest starbursts in our sample, each is V across. In the top left the white circle shows the size of the 
SDSS fibre through which the spectra are obtained. The white cross indicates the center of the image and approximate position of the fibre. 
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Figure A2. \' SDSS postage stamp images of the 20 oldest starbursts in our sample. 
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